A Note on the Diagonalization of the Discrete Fourier Transform 



o 



in 



C/2 



(N 



(N 
O 

o 



Zilong Wang*^'^ and Guang Gong^ 
^ School of Mathematical Sciences, Peking University, 
Beijing, 100871, P.R.CHINA 



^ ! ^ Department of Electrical and Computer Engineering, University of Waterloo 

Waterloo, Ontario N2L 3G1, CANADA 



Email: wzlmath@gmail.com ggong(Q)calliope.uwaterloo.ca 

June 5, 2009 



Abstract 



Following the approach developed by S. Gurevich and R. Hadani, an analytical formula of the 
canonical basis of the DFT is given for the case N = p where p is a prime number and p = 1 (mod 

V ■ 4). 

' Index Terms. Discrete Fourier transform, Weil representation, eigenvectors and orthonormal 

\^ ' basis. 
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1 Introduction 

The Discrete Fourier transform (DFT) has important appUcations in communication systems, and can 
be considered as an A'^-dimensional unitary operator F acting on the Hilbert space Ti = C(Zjv) by the 
formula 



- , — I 
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In the signal processing, the time domain and frequency domain are transformed by the DFT. A 
canonical basis, in other words, an orthonormal basis of eigenvectors for F will simplify the computation 
of the DFT. The main difficulty to get such a canonical basis is that F is an operator of order 4, and it 
has four distinct eigenvalues ±1, ±i with large multiplicity if the dimension > 4. The multiplicity of 
these eigenvalues depends on the value of n modulo 4, and was solved in P3|, although it was later shown 
to have been equivalent to a problem solved by Gauss in [S] . Unfortunately, no simple analytical formula 
for the eigenvectors is known. The research for finding different choices of eigenvectors, selected to satisfy 
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useful properties like orthogonality and to have simple forms, has been flourished in the literature |13| 
[S] H] U [ini, just listed a few here. 

A novel representation theoretic approach to the diagonalization problem of DFT in the case when 

= p is an odd prime number was introduced by Gurevich and Hadani in ]9^. This approach in- 
troduces the Weil representation [T7] of the finite symplectic group Sp = SL2{¥p) (will be precisely 
defined in Section 2) as the fundamental object of underlying harmonic analysis in the finite setting. 
More precisely, a centralizer subgroup of the DFT operator F in U (H.) (see definition in Section 2) is 
effectively described by using the Weil representation, which in this setting is a unitary representation 
p : SL2{¥p) U{'H) and the DFT is proportional to a single operator p{w) where w E SL2{Vp)- 
The centralizer subgroup of w is Tj„ which is a maximal algebraic torus (i.e., maximal commutative 
subgroup) in SL2{¥p). Then F commutes with p{Tio), and they share the same eigenvectors. 

By the above approach, a canonical basis (f>p of eigenvectors of the DFT and the transition matrix 
Qp from the standard basis to $p (discrete oscillator transform) for p = 1 (mod 4) were described by 
an algorithm in [3]. However, this algorithm has heavy computation cost, and the analytical formulas 
of eigenvectors of the Fourier matrix F were unknown. 

The vectors associated to the tori share many nice properties (see [5] [7] [H] for recent applications) 
and a simple analytical formula for the vectors associated to split tori was given in [16j . Based on [9] and 
|16j . in this paper, an analytical formula of the canonical basis of the DFT for the case of p = 1 (mod 
4) is given in Theorem 3, and their respective corresponding eigenvalues are determined in Theorem 
4. Then the discrete oscillator transform Op introduced in [5] can be obtained in a straightforward 
manner. 

The rest of the paper is organized as follows. In Section 2, we introduce the definitions of the one 
dimensional finite Heisenberg and Weil representations and the approach studying the eigenvectors of 
the DFT exhibited in [^. Then in Section 3, we give an analytical formula of the canonical basis of the 
DFT, and determine their respective corresponding eigenvalues. 

2 Preliminaries 

First, we introduce some basic concepts and notations which are frequently used in this paper. 

- For a given prime p, let 9 and rj denote the (p — l)th and pth primitive roots of unity in complex 
field respectively, i.e.. 



- We denote Fp as the finite field with p elements, and F* = Fp\{0} as the multiplicative group of 
Fp with a generator a. Then for every element b G F* there exist i with ^ z ^ p — 2, such that 
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b = a^. In other words, i = log^j b. 
- SL2{¥p) is the 2-diniensional special linear group over Fp consisting of all the 2x2 matrices 



- Let H = C(Fp) which is the p-dimcnsional Hilbert space containing all the function from Fp to C 
with the standard inner product, and U{Ti.) be the group of unitary operators on 7i. 

- A unitary representation of a group G on the Hilbert space 7i is a homomorphism p : G U (Ti.) 
which satisfies p{g ■ h) = p{g) ■ p{h) for Vg, h Q G. Specially, if G is an Abelian group, its 
representation p can be decomposed into a direct sum of 1-dimensional representation (character). 

2.1 The Heisenberg Representation 

Let {V,uj) be a two-dimensional symplectic vector space over the finite field Fp. For y{ti,Wi) £ V = 
Fp X Fp (i = 1, 2), the symplectic form uj is given by 



Considering V as an Abelian group, it admits a non-trivial central extension called the Heisenberg group 
H. The group H can be presented a.s H ^ V x ¥p with the multiplication given by 



It is easy to verify the center of i7 is Z = Z{H) = {g e H : gH = Hg] = {(0, 0, z) : z G Fp}. 

For a given non-trivial one dimensional representation (j) of the center Z, the Heisenberg group H 
admits a unique irreducible representation of H . 

Theorem 1 ( Stone- Von Neuman ) Up to isomorphism, there exists a unique irreducible unitary repre- 
sentation TT : _ff — > C/(7i) with central character (j), that is, n^z = 4' ■ 

The representation tt which appears in the above theorem is called the Heisenberg representation. In 
this paper, we take one dimensional representation of Z as 0((O, 0, z)) ~ rj^ . Then the unique irreducible 
unitary representation tt corresponding to (f> has the following formula 




such that ad ^ be ~ I where a, b,c,d £ F, 



{t2,W2)) = tiW2 - t2Wi. 



{ti,Wi,Zi) ■ {t2,W2,Z2) ^ {ti + t2,Wi W2, Zi Z2 + 2 ^Uj{{ti,Wi), (t2, ^2))). 



7r(t, w, z 



2 ^tW'\-z-\--wi 



(p{i + t) 



(1) 



for y^p e H, {t,w,z) £ H. 
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2.2 The Weil Representation 

The symplectic group Sp ~ Sp{V,ll}), which is isomorphic to SL2{Vp), acts by automorphism of H 

f a b \ 

through its action on the T^-coordinate, i.e., for V(i, w, z) G H and a matrix g = [ £ SL2{¥p), 

\c d J 

the action g on (t, w, z) is defined as 

g ■ {t, w, z) = {at + bw, ct + dw, z). (2) 

Let GL{TL) be the p— dimensional general linear group over C, and PGL{T-L) be the projective general 
linear group where PGL{Ti.) = GL{Ti) /C* . Due to Weil [l7j, a projective unitary representation 
p : SL2{¥p) PGL{T-C) is constructed as follows. Considering the Heisenberg representation tt : 
H U{Ti) and Vg G 5i2(IFp), a new representation is define as: tt^ : H ^ U{H) by tt^ (/i) = 7r(g(/i)). 
Because both tt and tt^ have the same central character (/), they are isomorphic by Theorem 1. By Schur's 
Lemma [14], TJom// (7r,7rf)^ C*, so there exist a projective representation p : SL2{¥p) PGL{Ti.). This 
projective representation p is characterized by the formula: 

pigMh)pig-') = ^{g{h)) (3) 

for every g G SL2{¥p) and h G H. A more delicate statement is that there exists a unique lifting of p 
into a unitary representation. 

Theorem 2 T/ie projective Weil representation uniquely lifts to a unitary representation 

p : SL2{¥p) U{n) 

that satisfies equation 

The existence of p follows from the fact [5] that any projective representation of SL2{Vp) can be 
lifted to an honest representation, while the uniqueness of p follows from the fact [9] that the group 
SL2{¥p) has no non-trivial characters when p 7^ 3. 

f a \ f 1 \ f 1 \ 

Note that S'L^iIP'b) can be generated by g„ = ] ■ Qt = \ 1 , and u> = 

^ ^ ^ ^ \ a-' J ' \ b I J' V -1 ; 

where a G F*, b G Fp. The formulae of their respective Weil representations for ga,gb and w are given 
in [7] as follows 

pigaMii) = a{a)ip{a-h) (4) 
p{gb)m) = v-'"''\i^) (5) 
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p — 1 

where (J : F* — > {±1} is the Legendre character, i.e., a{a) = a~2~ in Fp. 

Obviously, p{w) ~ F which is the DFT, and we denote p{ga) = Sa,p{gb) ~ Nf, for convenience. For 

V5 = [ " M eSL2{¥j,),iib^0, 
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c d y I (ad - 1)6-1 dylo fo-wlbd lM-1 ojl ab-^ 1 



Then the Weil representation of g is given by 

p{g)^SboNMoFoN,,-,. (7) 

If 6 = 0, then 
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a b \ _ a \ _ a OWlO 
c d \ c a-^ \ a-^ \ ac I 



Hence the Weil representation of g is as follows 

p{g) - Sa O Nac. (8) 
For more details about the Heisenberg and Weil representations, please see [B] [7] [IT] [T^ . 

2.3 Centralizer Subgroup of the DFT 

A. Maximal Algebraic Tori and T^, 

A maximal algebraic torus [2] in (IFp) is a maximal commutative subgroup which becomes 
diagonalizable over the original field or quadratic extension of the field. There are two classes of tori in 
SL2{Vp). The first class, called split tori, consists of those tori which are diagonalizable over Fp, while 
the second class, called non-split tori, consists of those tori which are not diagonalizable over Fp, but 
become diagonalizable over the quadratic extension Fp2 . 

Tm = \g : gw — wg,g G SL2(¥p)} is the centralizer group of w in SL2{¥p). It is easy to verify 

T„ = ||^^ ^ :a2 + 62^1,a,6eFp|. (9) 
If p = 1 (mod 4), then T^, is a split torus and conjugates to the standard diagonal torus 
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So is a cyclic subgroup of SL2{^p) with order p — 1. If p = 3 (mod 4), T^j is a non-split torus which 
is a cyclic subgroup of SL2{¥p) with order p + 1. 

B. Decomposition of Weil representation Associated with Tyj 

Because T^, is a cyclic group, restricting the Weil representation to T^,: p\T^^ : — > U {'H), we obtain 
a l-dimensional subrepresentation decomposition of p\T^ corresponding to an orthogonal decomposition 
of 7i(see [14] for basics of group representation theory). In other words, 

P\T^= X and (10) 

where At^^ is a collection of all the l-dimensional subrepresentation (character) x : — > C in the 
decomposition of Weil representation restricted to . 

If p = 1 (mod 4), X is the character given by x ■ ^p-i ^ C. We have dimTi.^ — 1 unless x — 
where a is the Legendre character of T, and dimTLa = 2. If p = 3 (mod 4), x is the character given by 
X ■ "Zp+i — > C. There is only one character which does not appear in the decomposition. For the other 
p characters x which appear in the decomposition, we have diniH^ = 1. 

Choosing a generator t £ T^, the character is generated by the eigenvalue x(t) of the linear operator 
p{t)^ and the character space naturally corresponds to the eigenspace of x(0- Because the eigenvalues 
of p{t) are almost different, it is easier to find a basis of orthogonal eigenvectors of p{t) than of the 
DFT. Since p{t) commutes with the DFT, the eigenvectors of p{t) are also the eigenvectors of the DFT. 
Thus, we obtain a canonical basis of the DFT. 

3 A Canonical Basis of the DFT 

If p = 1 (mod 4), Tw is a torus conjugating to the standard diagonal torus A = < : a e F* > , 

I V / J 

so which is a generator of T^, conjugates to = ( ) where a is a generator of F* i.e., there 

exist s € 5i2(Fp), such that t = sgaS~^ and p{t) = p{s)p{ga)p{s~^). Thus, the eigenvectors of p{t) can 
be determined by p{s) and the eigenvectors of p{ga)- In the following, we first present the results, then 
their proofs follow. 

Lemma 1 Let ga = i ) where a is a generator o/F* then {ip^ = {ipx{'i-)}o^i<p ■ ^ x < p} 

\ a^W 
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where 



1, for 2 = 



0, for i = 



■0o(«) = { and = < i for < x <p 

0, for i^O fl^'°g°', /or t 7^ 

Vp- 1 

is an orthonormal basis ofH and a collection of the eigenvectors of p{ga)- 

( 1 \ ,„ 1 _^ 

Lemma 2 Lei = I ^ j I where k = 4"; t/ien t = sffaS is a generator ofTyj. 

Thus $p = {ipx ■ = p{s)'ipx,0 ^ X < p} is a canonical basis of pit) = p{sgaS^^) and the DFT. 
More expUcitly, 



Theorem 3 Let 



^ r^^"^"'", forx = 



<Px{i) = < 



. P-i 

^ y 



Qxlog^JJ^a''U-^f-2 ^ aH^ ^ forO<X<p. 



Then $p — {(p^ ~ {V2;(*)}osSi<p : ^ x < p} is an orthonormal basis of Ti. and a collection of the 
eigenvectors of the DFT. 

Theorem 4 p^iO ^ x < p) is the eigenvector of the DFT corresponding to the eigenvalue {~i)^ where 
i = \J —\, i.e., 

Fipx = {-ifpx- 

Now we prove the above lemmas and theorems. Considering {5i : i £ ¥p} which is the orthonormal 
basis of Hilbert space H = C(Fp), where Si is defined as Si{j) = 5ij for Vi,j € ¥p, every vector 
ip = {p{i)} can be written as the form p = X^ieF Recall that SL2{¥p) can be generated by 

f a \ f 1 \ /0l\ 

Qa = I ] , Qb = \ and w = \ where a G F* and b e F„. Then their 

\0 a-^ J ^ \b 1 J \-l J " ' 

respective Weil representations (|4]),([5]), and ([6]) of pa, g^, and w can be rewritten as follows 

p{ga)5i = SaSi = a-{a)Sai (11) 
p{gb)S, ^ ^ 7j-^'"''" (12) 

piw)S, = FS, = ^ E (13) 
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Proof of Lemma 1. From pT]) . wc have 

P{ga)5t = Cr{,a)Sai = -Sat- 

Let Vi ^ V{Si),V2 = V{S2,S3r-- ,'5p-i), then it is obvious that H ^Vi^V2, < Vi,V2 >= 0, and 
p{ga){Vi) = Vi for i ~ 1,2. It is easy to see that p{ga)\Vi = —Id, so (5o is a eigenvector of p{ga) 
corresponding to the eigenvalue —1. The eigenfunction of p{ga)\V2 {x^^^ — 1), so the eigenvalues 
of p{ga)\V2 are e°,e\e^,--- ,6lP-2 which are different. We assert that 6I(t^-J' is the 

eigenvector associated to the eigenvalue 6^ {0 ^ j ^ p — 2), and it can be verified as follows 

p— 1 p— 1 



i=l 1=1 

i=l 

p-1 



Let X ~ — j. By normalizing the eigenvectors, we complete the proof. □ 
Proof of Lemma 2. Note that 



t ~ sgaS ^ 



1 2-^a'' \ ( a \ / -2-^a^ 



2-i(a-a-i) -2-ia''^(a-a-i) 
2-ia'^(a- a-i) 2-i(a-a-i) 



On the other hand, t conjugates to g^, so the order of < is p — 1. Thus, i is a generator of T^,. 

Proof of Theorem 3. Since t ~ sgaS^^, $p = {ipx = p{s)tpx : ^ a; < p} is a collection of the 
orthogonal eigenvectors of p{t) where (p^ and s are presented in Lemmas 1 and 2 respectively. 
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From (|12p . s has the following decomposition 



s ■ 



a'' 2-^ I \ 2a^'' / I 4~ia'= 1 / I -1 / I 2a^'' 1 



Then applying (HIl) , (Ull) , and for 1 x ^ p - 1, we have 



-1 fc\ P P-1 



T— ^EE^"'°'°''?"'-'"^''"'""'''^2-ia'=J (substitute I by 2-ia'=z) 





-1) 








-1) 


(t(2-i 





For 1' = 0, we have 

Because {tpx : ^ x < p} is an orthonormal basis and p(s) is a unitary matrix, {ipx : ^ x < p} 
is also an orthonormal basis of Ti.. Since p{s)F — Fp{s), ipx{0 ^ x < p) are not only the eigenvectors 
of p{s), but also the eigenvectors of the DFT. Note that a{2^^a'') is a constant, which completes the 
proof. □ 
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Proof of Theorem 4. It can be verified as follows, for a; = 0. we have 



1 

1 _ , P"^ 

= _^2 ia'=t^^^^2 la-'-jt^^ (substitute i - fl'^t by t) 



t=0 



1 k^2 

— if (Gauss sum) 



For a: 7^ 0, we have 

-ip-i 



p-i p-i 



-^=y e-'°^^-'"^y r^-''i-'"^--)"-^'"-''^"+-' (substitute 3 by a'^'j) 

n-kx P-1 P-1 



i=i t=o 



j=l t=0 
-^^^ j=l t=0 

|-^xiog,,,^a'^-o--t)^-2-ia''t= (Gauss sum) 



10 



Acknowledgment 

The authors would hke to thank Grcvich, Hadani and Sochen for their help during the course of 
conducting this work. 

References 

[1] N.M. Atakishiyev and K.B. Wolf. Fractional Fourier-Kravchuk transform. J. Opt. Soc. Am., Vol. 
14, No. 7, 1997, pp. 1467-1477. 

[2] A. Borel, Linear Algebraic Groups. Graduate Texts in Mathematics, vol. 126, Springer, New York, 
1991. 

[3] F.R. Beyl, The Schur multiplicator of SL{2, Z/mZ) and the congruence subgroup property. Math. 
Zeit, 191, 1986. 

[4] C. Candan, M. A. Kutay and H. M.Ozaktas, The discrete fractional Fourier transform, IEEE Trans. 
Signal Processing Vol. 48, No.5, 2000, pp. 1329-1337. 

[5] B.W. Dickinson and K. Steiglitz, Eigenvectors and functions of the discrete Fourier transform, IEEE 
Trans. Acoustics Speech and Signal Proc, Vol. 30, No. 1, 1982, pp. 25-31. 

[6] S. Gurevich, R. Hadani, and N. Sochen. The finite harmonic oscillator and its applications to 
sequences, communication and radar. IEEE Trans. Inform. Theory, Vol. 54, No. 9, September 2008, 
pp. 4239-4253. 

[7] S. Gurevich, R. Hadani, and N. Sochen, On some deterministic dictionaries supporting sparsity. 
Journal of Fourier Analysis and Applications, Vol. 14, No. 5-6, December 2008, pp. 859-876. 

[8] S. Gurevich, R. Hadani, and N. Sochen, Group representation design of digital signals and sequences, 
the Proceedings of the International Conference on Sequences and Their Applications (SETA), 2008, 
Sep. 14-18, 2008, Lexington, KY, USA. Sequences and Their Applications-SETA 2008, LNCS 5203, 
S.W. Golomb, et al. (Eds.), Springer, 2008, pp. 153-166. 

[9] S. Gurevich and R. Hadani, On the diagonalization of the discrete Fourier transform. Applied and 
Computational Harmonic Analysis, to appear 2009. 

[10] M.T. Hanna, N.P.A. Self, and W.A.E.M. Ahmed, Hermite-Gaussian-like eigenvectors of the discrete 
Fourier transform matrix based on the singular-value decomposition of its orthogonal projection 
matrices, IEEE Trans. Circ. Syst. I, Vol. 51, No. 11, 2004, pp. 2245-2254. 



11 



[11] S.D. Howard, A.R. Calderbank. and W. Moran, The finite Heisenberg-Weyl groups in radar and 
communications, EURASIP J. Appl. Signal Process, 2006, pp. 1-12. 

[12] R. Howe, Nice error bases, mutually unbiased bases, induced representations, the Heisenberg group 
and finite geometries, Indag. Math. (N.S.), Vol. 16, No. 3-4, 2005, pp. 553-583. 

[13] J. H. McClellan and T. W. Parks, Eigenvalues and eigenvectors of the discrete Fourier transforma- 
tion IEEE Trans. Audio Electroacoust, Vol. 20, No. 1, 1972, pp. 66-74. 

[14] J. P. Serre, Linear Representations of Finite Groups. Graduate Texts in Mathematics, Vol. 42, 
Springer, New York, 1977. 

[15] B. L. van der Waerden, Moderne Algebra, Springer, 1931. 

[16] Z. Wang, G. Gong, New sequences design from Weil representation with low two-dimensional 



correlation in both time and phase shifts, http:/ /arxiv.org/abs/0812.4487. Technical Report 2009-1, 
University of Waterloo, 2009. 

[17] A. Weil, Sur certains groupes d'operateurs unitaires. Acta. Math., Vol. Ill, 1964, pp. 143-211. 



12 



